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ABSTRACT 

^  A  software  package  for  processing  conductivity- 
temperature-depth  (CTD)  data  is  described.  The 
package  includes  features  for  editing,  correct¬ 
ing,  filtering  and  pressure  sorting  to  produce 
working  data  files  for  graphic  and  analysis  work. 
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A  COMPUTER  PROGRAM  FOR  PROCESSING 
CONDUCTIVITY -TEMPERATURE-DEPTH  (CTD)  DATA 


I .  INTRODUCTION 

CTD  (conductivity-temperature-depth)  data  acquired  with  a 
Neil  Brown  Instrument  Systems  CTD  system  (or  equivalent  data) 
require  certain  processing  steps  before  they  are  generally 
usable  for  oceanographic  analysis.  This  note  describes  a  proc¬ 
essing  software  package  (CTDPRG),  which  was  developed  at  the 
Naval  Oceanographic  Office  and  has  been  used  in  a  production 
mode  since  1978.  Some  minor  modifications  have  been  made  since 
that  time,  but  the  basic  structure  has  not  changed. 

Input  data  for  CTDPRG  must  be  in  the  form  of  a  Fast  Easy 
Binary  (FEB)  file  (Hallock,  1980)  and  must  consist  of  the 
physical  variables:  pressure  (dbar),  temperature  (°C)  and 

conductivity  (mmho/cm).  Hence,  a  pre-processing  step  is 
required  which  will  be  specific  to  the  original  data  recording 
format.  Output  is  quite  flexible  ranging  from  the  input  data 
itself  through  corrected,  pressure -sorted  and  subsampled  data 
with  calculated  quantities  such  as  salinity  and  density. 

Section  II  is  a  general,  functional  description  of  algo¬ 
rithms  used.  Section  III  provides  details  for  running  the 
programs  with  several  examples  of  run  setups.  Section  IV  con¬ 
tains  complete  program  listings. 

II.  GENERAL  DESCRIPTION 

CTDPRG  is  designed  to  handle  CTD  data  acquired  as  verti¬ 
cal  profiles,  horizontal  tows,  or  a  combination  of  these.  The 
program  was  developed  around  vertically  profiled  data;  how¬ 
ever,  and  much  of  its  complexity  goes  toward  converting  the 
time-base  CTD  data  to  functions  of  pressure,  when  this  conver¬ 
sion  is  required. 

The  following  list  is  a  summary  of  the  processing  steps 
performed  by  CTDPRG. 

•  Read  input  controls  from  cards  or  keyboard. 

•  Read  data  from  mass  storage  file. 

•  Correct  isolated  spikes  (wild  data  points). 

•  Apply  time-domain  filters  to  pressure,  temperature  and 
conductivity  (this  includes  T-C  matching). 

•  Correct  conductivity  values  based  on  water  sample 
comparisons. 

•  Apply  temperature-pressure  dependent  cell  correction. 

•  Separate  up-going  and  down-going  profiles. 

•  Pressure-sort  data  by  averaging  into  0.1  dbar  inter¬ 
vals. 


•  Low-pass  filter,  in  pressure,  and  subsample  to  mul¬ 
tiple  of  0.1  dbar. 

•  Compute  selected  quantities  (e.g.,  salinity,  sigma-t) 
from  P,  T,  C. 

•  Write  results  to  mass-storage  files. 

All  steps  (except  the  first  three  and  last)  are  optional 
and  can  be  selected  in  several  combinations. 

Input  Data  Sets 

The  basic  unit  of  processing  is  the  input  data  set.  This 
is,  typically,  a  single  CTD  cast  (down  and  up,  if  both  were 
recorded).  It  might  be,  however,  a  series  of  down  and  up 
profiles  acquired  during  one  instrument  deployment.  Usually  an 
input  data  set  corresponds  to  a  single  input  file,  but  this  is 
not  required.  There  may  be  more  than  one  data  set  per  file  or 
a  data  set  may  span  several  files.  All  this  is  handled  by  the 
input  controls.  The  input  data  set  is  in  a  FEB-file  containing 
time,  pressure  (dbar),  temperature  (°C)  and  conductivity 
(mraho/cm).*  Additional  variables  could  be  carried  with  minor 
modifications.  The  last  FEB-segment  in  the  data  set  must  have 
ID0CR(1)=1;  preceding  segments  (for  multi-segment  data  sets) 
must  have  IDOCR(1)=0.  Input  files  are  thus  prepared  by  a 
pre-processing  program  that  is  tailored  for  the  particular 
instrument  and  deployment  mode. 

Time-Domain  Filters 

One  of  the  most  significant  problems  with  CTD  data  has 
been  the  effects  of  the  response  mismatch  between  temperature 
and  conductivity  sensors.  A  number  of  attempts  have  been  made 
to  correct  this  through  the  design  and  application  of  filters. 
The  most  widely  used  method  is  that  described  by  Fofonoff  et 
al.  (1974),  which  assumes  an  exponential  response  for  the 
temperature  sensor.  More  recently,  as  a  result  of  work  done  at 
the  Naval  Oceanographic  Office  (N00)  involving  direct  measure¬ 
ment  of  instrument  response,  specific  filters  have  been 
designed  for  temperature  and  conductivity,  based  on  experimen¬ 
tally  determined  sensor  response  functions  (see  Mayoral, 
1982).  In  CTDPRG  there  are  three  options:  no  T-C  filtering, 
filtering,  the  three-point  Fofonoff  method  (the  user  specifies 
only  the  time-constant) ,  or  explicitly  supplied  filters  of  the 
user's  choice.  These  are  applied  to  the  raw  data,  in  the  time 
domain,  prior  to  other  processing. 

A  smoothing  filter  is  automatically  applied  to  pressure 
data  to  remove  high-frequency  sensor  noise. 

Conductivity  Calibration 

With  the  newer  conductivity  cells  on  the  MK  III  CTD 
systems,  coupled  with  proper  cleaning  procedures  and  frequent 
laboratory  calibration  adjustments,  conductivity  data  seldom 


♦See  Hallock  (1980)  for  details  on  FEB-file  structure. 
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has  to  be  adjusted  from  the  recorded  values.  However,  if 
salinity  checks  derived  from  water  samples  indicate  an  offset 
from  proper  calibration,  a  multiplicative  correction  can  be 
supplied  by  the  user.  There  is  also  a  manufacturer-supplied 
correction  for  pressure  and  temperature  distortion  which  is 
applied  following  the  filtering  step. 

Upcasts  and  Downcasts 

If  no  editing  based  on  pressure  is  to  be  done  (i.e.,  no 
identification  of  upcasts,  downcasts  or  constant  level  data), 
following  conductivity  correction,  selected  additional 
variables  are  computed,  and  an  output  data  set  is  written  that 
contains  the  same  number  of  scans  (cycles)  as  the  input. 


A  second  option  is  to  retain  only  upcasts  and/or 
downcasts,  based  on  user  defined  pressure  limits,  which  are 
then  written  as  separate  output  data  sets,  following 
computation  of  additional  variables. 

A  third  option,  which  necessitates  the  second,  is  to 
convert  each  upcast/downcast  to  a  pressure  series  as  opposed 
to  a  time  series.  This  case  assumes  that  the  instrument 
descent /ascent  rate  is  2  m/sec  or  less.  (If  this  value  is 
exceeded  in  the  mean,  some  software  modifications  would  be 
necessary.)  At  a  nominal  descent/ascent  rate  of  1  m/sec,  about 
three  scans  of  data  are  acquired  every  0.1  dB.  If  pressure 
sorting  is  selected,  time,  temperature  and  conductivity  values 
in  0.1  dbar  intervals  are  averaged,  creating  a  pressure  series 
at  that  resolution.  Retrograde  (in  pressure)  data  are  re¬ 
jected.  The  pressure  series  is  then  subsampled  at  some 
user-defined,  integral  multiple  of  0.1  dbar  following  the  ap¬ 
plication  of  a  low-pass  filter  to  prevent  aliasing.  The  filter 
is  automatically  generated,  and  its  cutoff  wavenumber  is  the 
new  Nyquist  wavenumber  of  the  subsampled  series.  At  this  point 
additional  variables  are  computed  and  data  are  written  to  out¬ 
put  files.  The  third  option  allows  the  user  to  select  a  degree 
of  smoothing/subsampling  commensurate  with  his  requirements. 

Computed  Variables 

Including  the  four  input  variables,  up  to  twelve  are 
available  as  output.  The  user  may  select  any  subset  of  these. 
The  eight  computed  variables  are  as  follows: 

•  Salinity  (parts  per  thousand  [ppt])  or  S:  computed  with 
pressure,  temperature  and  conductivity  ratio,  according  to 
Bennett  (1976),  Dauphine  and  Klein  (1977)  or  Lewis  (1980). 

•  Sigma  or  «  :  the  in  situ  density  of  a  water  parcel.  Com¬ 
puted  with  pressure,  temperature,  and  salinity  according  to 
Fofonof f  (1962).  (Sigma  =  (density  (g/cc)™l . )*1000.  ) 

•  Sigraa-t  or  the  density  of  a  water  parcel  with  in 
situ  temperature  and  salinity  but  at  atmospheric  pressure. 

•  Brunt-Vaisala  Frequency  (cycles  per  hour)  or  N:  computed 
with  pressure,  temperature,  and  salinity  according  to  Fofonoff 


(1962).  Essentially  a  local  potential  density  derivative.  The 
convention  is  adopted:  N  ■  V]NZ |  sgn  (N^)  to  allow  for  in¬ 
stabilities. 

•  Sound  Speed  (m/sec):  computed  with  pressure,  tempera¬ 
ture,  and  salinity  according  to  Wilson  (1960). 

•  True  Depth:  for  most  purposes,  pressure  in  decibars  can 
be  interpreted  as  depth  in  meters;  however,  this  assumes  an 
average  water  density  of  1.02  6/cc.  If  depth  accurate  to 
more  than  a  few  percent  are  required,  this  variable  can  be 
selected.  It  is  an  approximation  to: 

-  -  /  It  dp 
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where  a  -  specific  volume, 

g  =  acceleration  of  gravity, 
p  =  pressure, 
z  ■  depth  below  surface, 

•  Potential  Temperature  (°C)  or  9:  The  temperature  that 
a  parcel  of  water  would  have  if  brought  adiabatically  to  the 
surface.  Computed  with  pressure,  temperature  and  salinity 
according  to  Fofonoff  (1962). 

•  Potential  Density  or  *q:  using  0  rather  than 

in  situ  temperature;  the  density  of  a  water  parcel  brought 
adiabatically  to  the  surface.  (  <rg  is  calculated  using  the 

routine  with  9  from  the  above  routine.) 

A  typical  standard  processing  product  might  be:  retain 
only  downcasts  at  1  meter  resolution,  and  store  pressure,  tem¬ 
perature,  salinity  and  <ff  This  provides  a  compact  version 
of  data  for  hydrographic  analysis  work.  The  best  "standard" 
file  is  likely  to  be  different  for  each  user  group;  thus,  it 
is  recommended  that  the  input  files  (edited,  raw  data)  be 
considered  the  primary,  archived  data  base. 

III.  DETAILS  OF  OPERATION 

This  section  consists  of  two  parts:  a  tour  of  CTDPRG  and 
its  subroutines;  a  tutorial  on  setting  up  input  and  running 
the  program. 

Functional  Program  Description 

Figure  1  is  a  block  diagram  of  all  the  routines  in  the 
package.  CTDPRG  is  the  driver,  whose  functions  are  to  read  all 
input  controls,  to  set  up  filters  to  be  used  on  input  data, 
and  to  initiate  the  processing  of  each  input  data  set  by  call¬ 
ing  ACCESS. 

GENWHF  generates  a  three-point  deconvolution  filter  for 
temperature  based  on  a  user-specified,  sensor  time-constant. 
This  is  described  in  detail  by  Fofonoff  et  al.  (1974). 


GENFLT  is  a  general  filter  generating  routine  (see 
Brooks,  1976).  It  is  called  in  CTDPRG  to  produce  a  low-pass 
filter  for  pressure  data  to  remove  sampling  noise. 

ACCESS  controls  the  majority  of  the  processing.  Input 
data  cycles  (scans)  are  accessed,  one  at  a  time,  via  calls  to 
GETREC.  If  pressure  sorting  is  not  selected,  output  segments 
are  prepared  using  a  user-specified  combination  of  variables 
from  the  input  as  well  as  quantities  computed  via  DERIVE. 

When  upcast/downcast  separation  is  specified,  each  up¬ 
cast/downcast  becomes  a  separate  data  set  and  is  written  to 
unit  3/2.  Otherwise,  all  output  goes  to  unit  2  as  one  data 
set. 

With  pressure  sorting,  which  necessitates  upcast /downcast 
separation,  input  data  cycles  in  0.1  dB  intervals  are  aver¬ 
aged,  producing  a  pressure  series  at  that  resolution.  Retro¬ 
grade  data,  which  result  from  ship  motion  or  winch  reversals, 
are  rejected.  For  a  nominal  descent  rate  of  1  meter  per  second 
there  should  be  about  3  cycles  per  interval.  If  an  interval 
should  turn  up  empty,  it  is  filled  with  data  from  the  previous 
interval.  If  descent  rates  systematically  exceed  1.5-2  m/s 
this  interval  should  be  increased  accordingly. 

The  pressure  sorted  data  can  be  subsampled  at  integral 
multiples,  say  n,  of  0.1  dB.  When  this  is  done,  ACCESS  calls 
GENFLT  to  produce  a  low-pass  filter  which  is  applied  to  the 
0.1  dB  data,  centered  at  every  nth  cycle.  The  cutoff  of  the 
filter  is  at  the  new  Nyquist  of  the  subsampled  series  and  pre¬ 
vents  aliasing.  For  example:  n  =  10  implies  subsampling  every 
1  dB;  cutoff  (half-amplitude  point)  is  at  0.5  cycles/dB;  the 
filter  is  applied  only  where  samples  are  to  be  extracted 
(e.g.,  every  10th  point),  thus  avoiding  unnecessary  computa¬ 
tion. 

Following  pressure  sorting  and  subsarapling,  DERIVE  is 
called  for  selected  computed  variables.  Segments  are  then 
assembled  and  written  to  the  output  file(s). 

Subroutine  GETREC  is  the  input  program.  It  reads  segments 
from  the  input  file  (unit  4),  applies  conductivity  cell  corr¬ 
ections,  pressure  smoothing  filters,  conductivity-temperature 
matching  filters  and  pressure  limits.  Data  are  then  passed  to 
ACCESS,  one  cycle  at  a  time.  Upper  and  lower  pressure  limits, 
specified  by  the  user,  are  employed  by  GETREC  to  distinguish 
upcast  and  downcast  d^ta. 

NBCCC  applies  a  manufacturer  supplied  cell  correction  to 
conductivity  which  compensates  for  pressure  and  temperature 
distortions  of  cell  geometry.  It  also  applies  a  user-specified 
multiplicative  conductivity  adjustment  (usually  =  1.)  which 
corrects  for  salinity  offsets  determined  by  field  comparisons. 

Time-domain  filters  are  described  in  Section  II. 
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When  an  end-of-data  set  (IDOCR(l)  ■  1)  is  detected,  con¬ 
trol  returns  to  the  main  program  and  new  input  cards  are 
sought.  If  an  end-of-file  is  encountered,  however,  GETREC 
calls  IUSET  which  attempts  to  dynamically  attach  another  user- 
specified  file  name  to  unit  4.  If  no  file  name  is  supplied, 
the  program  terminates.  (On  some  machines  dynamic  control 
commands  may  not  be  permitted,  thus  requiring  one  execution 
per  input  file.) 

Subroutine  DERIVE  accepts  a  cycle  of  pressure,  tempera¬ 
ture,  conductivity  and  time  from  access  and,  after  calling  the 
appropriate  variable  computation  routines,  restructures  the 
cycle  according  to  user  specifications.  The  variable  computa¬ 
tion  routines  are  described  in  Section  II. 


Subroutine  ZREAD  and  ZWRIT  (described  by  Hallock,  1980) 
are  input  and  output  routines  for  FEB-files.  A  call  to  ZREAD/ 
ZWRIT  transfers  one  segment  of  data  between  main  memory  and  a 
mass  storage  file.  Typically,  a  data  set  (e.g.,  CTO  profile) 
consists  of  several  FEB  segments;  each  segment  is  made  up  of 
several  header  blocks  followed  by  a  data  block.  Segments  are 
written  sequentially  to  a  mass-storage  (disk)  file  by  nonfor- 
matted,  FORTRAN  I/O  statements. 


Input  Control  Cards 

Essential  information  for  setting  up  and  running  CTDPRG  is 
provided  through  internal  documentation  (comment  cards),  which 
are  presented  below. 


Input  Card  1:  Printout  identification — up  to  48  Alphanu¬ 
meric  characters.  This  can  be  anything. 


Input  Card  2:  N,  LN,  NWMAX,  KCHOP,  NINFIL  (Free  Format). 

N:  Subsampling  interval  for  pressure  sort¬ 
ed  data.  Integer  multiples  of  0.1  dB, 
e.g.,  N=20  results  in  data  every  2  dB. 
N=0  results  in  no  pressure  sorting. 

LN:  Sharpness  factor  for  subsampling  fil¬ 
ter.  This  is  an  integer  equal  to  1  or 
greater.  Larger  values  require  more 
computer  time,  (ra  =  LN*N  =  half-width 
of  filter.  See,  GENER1  in  Brooks,  1976, 
for  a  more  detailed  explanation.) 


NWMAX:  Maximum  no.  of  cycles  in  output  seg¬ 
ments.  Typically,  1000. 


KCH0P=1:  Process  all  data.  No  pressure  sorting 

or  subsampling.  Pressure  limits  are  not 
applied.  Output  goes  to  unit  2. 

KCH0P*2:  Pressure  sort  and  subsaraple  by  N  (un¬ 
less  N  ■  0).  Output  downcasts  to  unit 
2,  upcasts  to  unit  3. 


KCH0P*3:  Same  as  2  but  output  only  downcasts. 
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KCH0P=4:  Same  as  2  but  output  only  upcasts. 

NINFIL:  No.  of  input  file  names  to  be  supplied. 
If  =0,  one  input  file  is  assumed  on 
unit  4.  Otherwise,  supplied  names  are 
dynamically  attached  to  unit  4  by  IU3ET 
(presently  configured  for  UNIVAC). 

Input  Card  3:  (IND(I),  1=1,  12)  (Free  format)  output  var¬ 
iable  selection.  A  string  of  12  integers  (1 
for  yes,  0  for  no)  indicating  variables  to 
be  output  from  the  following  list:  Time, 
pressure,  temperature,  salinity,  sigma-t, 
sigma,  Brunt-Vaisala  frequency,  sound 
speed,  depth,  potential  temperature,  poten¬ 
tial  density. 

Input  Filename  Cards:  up  to  24  alphanumeric  characters, 

left-justified.  If  NINFIL=0,  no 
cards  are  expected. 

The  following  group  of  cards,  the  data  set  card  and  pres¬ 
sure  limit  cards,  is  repeated  for  each  input  data  set  to  be 
processed.  . 


Input  Data  Set  Card:  NMSER,  NPLIMS,  CFTR ,  TIMCON,  ISAL 

(Free  Format) 

NMSER:  Sequence  no.  of  input  data 

set  to  be  processed.  (Usual¬ 
ly  this  is  equal  to  1,  but  if 
there  were,  say,  n  data  sets 
to  be  skipped  beginning  in 
the  first  input  file,  NMSER 
would  be  set  to  n+1.) 

NPLIMS:  No.  of  pressure  limit  cards 
for  this  data  set. 


CFTR:  Conductivity  cell  correc¬ 
tion  factor  (default=l). 

TIMCON>0:  Temperature  sensor  time  con¬ 
stant  divided  by  sampling 
time  interval.  WHOI  (Fofo- 
noff  et  al.  temperature 
deconvolution. ) 

=0:  No.  temperature  or  conductiv¬ 
ity  filtering. 

*-l:  Use  same  filters  that  were 
used  for  previous  data  set. 
(Not  permissible  for  first 
data  set . ) 


=-5:  Read  special  filters  from 
cards  (unit  5)  following 
pressure  limit  cards. 
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All  Else:  Read  special  filters  from  file 
on  unit  8. 

ISAL  =  0:  Salinity  computed  according  to 
Lewis  and  Perkin  (1980) 

=  1:  Salinity  computed  according  to 
Bennet-Dauphinee. 

Pressure/Cycle  Limit  Cards:  PTOP,  PBOT,  DELTAP,  II,  12 

(Free  format) 

PTOP,  PBOT  are  upper  lower  pressure  limits  to  be  applied 
to  input  data  cycles  between  cycles  II  and  12.  Each  sequence 
of  cycles  spanning  this  pressure  interval  is  made  into  an 
output  data  set.  Pressure  limits  are  irrelevant  for  KCH0P=1. 

DELTAP  is  a  slop  factor  (usually  about  2  dB)  that  ’  ’  em¬ 
ployed  to  prevent  small  pressure  fluctuations  from  trip*  t  up 
the  pressure  limit  logic.  It  ensures,  for  example,  it  a 
downcast  has  ended  before  the  subsequent  upcast  is  rt  ;ni- 
zed.  For  single-cast  data,  DELTAP  can  be  set  to  zero. 

II,  12 — There  is  a  global  cycle  index,  IE,  which  is  to 

1  at  the  beginning  of  an  input  data  set.  It  is  increment  -  -or 
each  cycle  encountered  within  the  data  set.  (An  input  daLa  set 
typically  spans  more  than  one  FEB  segment  and  occasionally 
more  than  one  file.)  The  first  profile  start,  defined  by  PTOP, 
is  not  sought  until  IE  reaches  II.  After  each  profile  is  com¬ 
pleted,  if  IE  has  reached  12,  the  next  pressure  limit  card  is 
invoked.  When  the  cards  are  exhausted  or  the  input  data  set 
ends,  the  next  data  set  card  is  read.  This  feature  allows 
leading  and  trailing  data  within  a  data  set  to  be  skipped, 
which  can  save  computer  time  or  reject  defective  data. 

Usually,  only  one  pressure  limit  card  is  required  for  each 
input  data  set,  particularly  for  a  single  profile.  However, 
for  a  series  of  repeated  profMes  which  are  included  in  a 
single  data  set  (i.e.,  a  "yo-yo"  cast),  it  may  be  necessary  to 
change  pressure  limits  several  times  to  maximize  data  recovery 
and  to  reject  sections  of  bad  data. 

In  the  simplest  case,  for  a  single  cast  data  set,  II  can 
be  set  to  1  and  12  to  a  large  number  which  exceeds  the  number 
of  cycles  in  the  data  set,  and  everything  will  be  processed. 
(Default  values  of  1  and  1,000,000  result  for  II  =  12  =  0.) 

Special  Filter  Input:  NWTT ,  NWTC,  LAGT ,  (WT(I),  1=1, 

NWTT ) ,  WC(I),  I=NWTC)  (Free  Format) 

These  are  read  from  cards  if  TIMCON  =  -5.  If  TIMECON  =  -2, 

-3,  -4,  -6,  -7,...  they  are  read  from  a  file  on  unit  8.  If 

TIMCON  >  -2,  this  input  is  not  expected.  (See  TIMCON 

above. ) 

NWTT,  NWTC:  No.  of  temperature,  conductivity,  filter 
weights . 
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LAGT:  No.  of  cycles  to  lag  temperature  relative  to  pres¬ 
sure  and  conductivity.  (The  lag  might  be  built  into  fil¬ 
ter,  in  which  case  LAGT  =  0.) 

WT,  WC:  Temperature,  conductivity  filter  weights. 

Examples  of  Input  Controls 

Example  1:  Three  input  files,  each  containing  a  downcast 
and  an  upcast  between  about  8  m  and  2000  m  depth.  Only 
downcasts  are  to  be  processed.  WHOI  (Fofonof,  et  al.,1974) 
temperature  deconvolution  is  to  be  used — time  constant  is 
100  ms,  sampling  interval  is  32  ms.  Output  pressure  reso¬ 
lution  is  to  be  1  dbar.  Output  variables  are  to  be:  P,  T, 
S,  0% ,  N.  Salinity  to  be  computed  by  practical  salinity 
scale  '78. 


Card  1 
Card  2 
Card  3 
Card  4 
Card  5 
Card  6 
Card  7 
Card  8 


THIS  IS  EXAMPLE  1. 
10,1,1000,3,3 
0,1, 1,0,1, 1,0, 1,0, 0,0,0 
FILE1 
FILE2 
FILE3 

1,1,1,3.125,0 

10,2000,1,1,0 


Cards  7  and  8  are  repeated  for  Files  2  and  3. 

Example  2:  One  input  file  containing  a  continuous  record 
of  two  downcasts  and  two  upcasts  (i.e.,  one  instrument 
deployment).  The  instrument  was  ir  the  water  for  one  hour 
before  the  first  downcast  began  and  again  following  the 
final  upcast.  Thus,  the  first  usable  data  scan  is  about 
no.  112,000.  The  four  profiles  took  about  one  hour  so  the 
last  usable  scan  was  around  250,000.  Temperature  deconvo¬ 
lution  filters  are  to  be  read  from  unit  8;  a  1.01  conduc- 
is  to  be  made;  salinity  is  to  be  computed  by  Bennet- 
Dauphinee  algorithms.  Data  are  to  be  subsampled  at  50  cm. 
Both  upcasts  and  downcasts  are  to  be  saved,  between  100 
and  400  dbar. 


Card  1:  THIS  IS  EXAMPLE  2. 

Card  2:  5,1,1000,2,0  (input  assumed  on  unit  4) 

Card  3:  0,1, 1,0, 1,1,0, 1,0, 0,0,0 

Card  4:  1 , 1 , 1 .01 , -3 , 1 

Card  5:  100,400,2,112000,250000 

Example  3:  One  input  file  containing  500  downcasts  and  up¬ 
casts  between  about  90  and  210  dbar  ("yo-yo"  cast).  Data 
are  to  be  processed  as  in  example  1,  saving  downcasts  at 
1  dbar  resolution.  Prior  examination  of  the  data  reveals 
that  downcasts  201-208  (scans  1.2  x  10®  -  1.248  x  10®) 
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extend  only  to  190  dbar:  three  sets  of  pressure  limits  are 
thus  required,  to  recover  the  maximum  amount  of  data. 

Card  1:  THIS  IS  EXAMPLE  3. 

Card  2:  10,1,1000,3,0 

Card  3:  0, 1 , 1 , 0 , 1 , 1 , 0 , 1 , 0 , 0 , 0 , 0 

Card  4:  1,3,1,3.125,0 

Card  5:  95,205,2,1,1200000 

Card  6:  95,185,2,1200000,1248000 

Card  7:  95,205,2,1248000,4000000 

Example  4:  Input  file  consists  of  three  downcast-upcast 
pairs.  For  diagnostic  purposes,  it  is  necessary  to  process 
downcast  data  at  original  resolution  (i.e.,  no  pressure 
sorting).  Output  variables  are  to  be:  time,  pressure,  tem¬ 
perature,  conductivity  and  salinity.  Pressure  limits  are 
10  to  500  dbar.  Special  filters  for  temperature  and  con¬ 
ductivity  are  to  be  supplied  on  cards. 

Card  1:  THIS  IS  EXAMPLE  4. 

Card  2:  0,0,1000,3,0 

Card  3:  1,1, 1,1, 1,0, 0,0, 0,0, 0,0 

Card  4:  1,1,1, -5,0 

Card  5:  10,500,0,1,0 

Card  6:  5 , 5 , 0 , . 1 , . 2 , .4 , .2 , . 1 , .0625 , . 125 , .625 

Card  7:  .125, .0625 

These  are  representative  examples  of  how  CTDPRG  may  be 
used.  Many  other  combinations  are  possible.  Some  experimenta¬ 
tion  may  be  necessary  for  unusual  requirements. 
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CONTROLS 


Functional  Block  Diagram  of  CTD  Processing  Programs 


APPENDIX 
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t.CWHt 

ccM»(iieiu:i> 


c  risk  era  ntctitmi  ntutn 

c 

C  INPUT  F1UI  (UNIT  *1  Alt  PROPERLY  FI(MHQ  FES  FILES  OF 

C  IK  CTO  on*.  «HIilU  AC  ((Ml  PRESSURE,  rC*K»iUl*tt 

t  CONOUCf I«IIY.  NNFI,  MH  SHOULD  »C  COUlSt  *MD  SILI10N 

C  NUMRII,  RESPECTIYELT,  UtUlCi  NARKS  INC  (MB  OF  *M 

C  IHFU  T  0*U  ill  1CAST1.  f«U  CAN  RE  HBOC  1H»N  ONt  IK  PUT 

c  rut  pci  on  *  set  tstt  celoui.  output  oocs  to  units 

C  2  A  HO  3  «  UNScPLIITCO  0*1*  <MCMOP=|f  4*0  D(**CAj7S  SO  >0  UNIT 

C  2  NH1LC  UPCASTS  CO  TO  UNIT  S.  TOO  FURTHER  DETAILS  SCC 

c  iitcinIl  DocuHCnmoN.  Input  conhol  cards  arc  oes- 

C  CSIKO  It  LOU. 

c 

CCPNON  FUN*  I/N1N0  , 1  NO  1201  t*CH9P  «NPL  1HS,»L  11,1001  .1PL1HI2,  TOO*  » 

•  KNHAX, JPXSCTI207, ICOLU.NMSER 
COMMON /SAL 1NC/*S*L 
CONNON/IOFlLC/falNFlL.lFCBIial 

CCMMON/CALlS/MP,  NT  ,NC ,t»C1,  NPt  *01  ,UT|*OI,UCt  *0  I  .POIF. 

•  IDTF  .CO  IF.Cf  TO 

ell C NS  ION  L*R£LI20l,lRUNItl 

C*T*  ILAICUII.IS1.12I  /’11ME‘, ’PRESS’, *TCNP’, ’COHO’, ‘SAl‘»‘S1ST’, 

•  ’Sl&STP  *,*S»FCPM’ .’SHOVEL*, ‘TOEPTM’, ’TCNPOT ’.’SISPOT’/ 

CAT*  NP,U«,JC.JC.F/10.2.t.l,.l/ 

CAT*  ISM/I/ 

0*1*  P0IF,10IF,C01F/.S,.2,.2/ 

C 

C  •«•*«*  INPUT  C*«0  Is  PltNTOUT  IDENTIFICATION.  UP  TO  *0  CHAOS. 

C  CAN  SC  A  IN  THINS  -  EVEN  CLANK  CADD. 

C 

ICAOIS. 1020, C NO SCOO 1  INUN 
1020  FCINaT 11**1 

KRlll  IN ,1021  llOUM 

10)1  fCRMATI IHl $*  CTO  POOCCSSlNS  OUNt  ’,***//> 


C  •«••*•  INPUT  CARD  2*  N, IN, NUN  II  ,*CHOP  ,N|NF  IL 
C 

C  »<  SUSSAHPLlNG  INTERVAL  FOR  PRESSURE  SORT  INS. 

C  INTEGER  WLUPLCS  OF  .1  31.  I.E.  N-20  RESULTS 

C  IN  OAT*  EYERY  2  OS.  ISO  RESULTS  IN  NO  SOUTHS. 

C 

C  LN,  SHARPNES  FACTOR  FOR  DCCINATION  FILTER. 

C  THIS  IS  HI  INTEGER  EQUAL  TO  1  OR  GREATER. 

C  IT  OEURhInCs  THE  NO.  OF  HEIGHTS  *N0  THUS 

C  INC  SHARPNESS  OF  ThC  DCCINATION  FILTER. 

C  NNNAAs  RAX  HO.  V  CYCLES  JN  OUTPUT  SEGMENTS, 

t  IAS  A  RULE,  THIS  SHOULD  SE  ASOUT  1000. 1 

C  KCHOP- 1 s  PROCESS  ALL  OATAj  NO  PRESSURE  sorting 

c  OR  SURSanPLInC.  PRESSURE  LIMITS  ARC 

C  NOT  APPLIED, 

t 

t  «CM0P:2t  OUTPUT  UPCASTS  AND  OONNCASTSJ 

C  PRESSURE  SORT  AND  $UIS|N*l;  RT  N. 

C 

C  KCHOP* 2>  SANE  AS  2  BUT  OUTPUT  ONLY  DOUNCASTS. 

C 

C  KCHOP  I* »  SANE  AS  2  CUT  PROCESS  ONLY  UPCASTS. 

C 

C  RINFIU  NO,  OF  1»0T  FILENAMES  TO  SE  SUPPLIEO. 

C  IF  MINFILSO,  ONE  INPUT  FILE  IS  ASSUMED  ON 

C  UNIT  «.  OTHERNISC.  SUPPLIEO  NAMES  ARC  ASS1SNE0 

C  NYMAM1CAU.Y  TO  UNIT  «  SV  SR  IUSET.  THIS  NAY 

C  NOT  CORK  ON  Ul  MACHINES. 

C 

StAOIS,IOOO>ENO:«Y«IN,LN,NyNAI,KCH*P,NlNFIL 
IOCO  FORNATTI 

klllt  u,loqqtn,lm,nnmax,kch«p 
KCHOP  I  "KCHOP 

LOCO  FcRmATIIX,’N*’,1S,1X •  *LNS‘, IS ,  ll.’NuNAXs’  , 

cis,ik,’nchopx*,isi 

c 

C  •  INPUT  CARO  s*  OUTPUT  YUIIaSLC  SELECTION.  I  STRING  OF 

C  12  INTEGERS  II  FOR  YES,  0  FOR  NOI 

C  INDICATING  YA RUBLES  TO  RE  0UTPU1,  FROM 

C  TM(  fOLLGUpC  LIST, 

c 


C  TINE. PRESSURE, lENPERATURCfCONOUCTlVITV, 

c  uiMinitiMi'ii  iimit  m  netucwi, 

c  mimo  mto.  true  oepin,  potential 

C  1dM«  M(«rOR<UU  MMI1I. 

c 

c 

»:S 

BO  II  ls| ,20 

lFllNDIl».£0.g|00  tO  || 

4-J  *  1 

IMidIJKllllLUl 
II  CONTINUE 
AINOIJ 
C 

biIIE|A,*OOS» 

ioci  formati//1s, ‘variables  requested*! 

NRIICIA,AOOAIIIPMSET  IRI«RA|  <JI 
0000  rORMTIIlI.IQUO.lIIII 

•  «•••«  INPUT  FILENAME  CIUI  irOMlT  0041 

one  c»ro  rm  each  input  filename  or  up  to 
2*  characters  iincluoinb  •».  ir  ninfil=s, 
no  eons  are  tvccns. 

CO  Of]  Isi.o 
Oil  IFESIII  :•  • 

IFtOlMFlL  .CO. 01  00  TO  000 
CO  t«l  I=l,NlNFIL 
I|:|I-II***1 
12:11*1 

VSl  SCA0|S,100!,EN0=V9Vl  IlFEtl  nl  ,II=U«I2> 

loot  nniutfii 

9 SO  CONTINUE 


THE  FOLLOHlNS  BROUP  OF  ClHl)  THE  DOTO  ill 
"CARO  I  NO  PKSSURE  LIMIT  CMOt,  10  REPEATED 
FOO  EACH  llPNT  OOTO  OET  It  BE  PROCESSED. 

•  INPUT  OOTO  SET  COROl  NMSER  .NPLIMS.CFTR.T  INCON  .ISOL 


OH  SCR 

SEQUENCE  MO.  OF  INPUT  OOTO  SET  TO  BE  PNOCCSSE0. 

OPLIRS 

NO.  OF  PRESSURE  limit  CORDS  for  THIS  OOTO  SET. 

CF  IN 

CONOUCTIVITV  CELL  FACTOR  IDCFOULTsl.l 

IIpCOn 

ST  01 

TEMPERATURE  TINE  CONSTONT/SONPLINB  tine  interval 
UNO!  TEMPERATURE  DECONVOLUTION. 

S  Ol 

NO  FILTERS  OR  TEMP  OR  CORO. 

:  -II 

USE  FILTERS  FOR  PREVIOUS  DATA  SET.  INOT 
PERMISSIBLE  ON  FIRST  DATA  SCT.l 

s  -S« 

READ  SPECIM.  FILTERS  FROM  CAROS  FOLLOUINS 
PRESSURE  LIMIT  CORDIS*. 

ELSE! 

ISOL 

:  Ol 

’  salinity  computed  Ian  pssto. 

:  It 

salinitv  computed  Ion  rennet -oauphinee. 

IS  AtA0IS,1000,EN0:«VVINNSER,NPLINS,CFTR,TINC0N,lSAL 
NSOL zt 

IF  <ISOL.NE.O|MSOt:| 

PHI 
PC  si 
OUTTsI 
OhlCI 
NT IIIS1. 

NCIIISI. 

IFILH2 

IF  IIIRcBN.LT. 0.1  IFILTSB 
IFttlNCON.LT.-l.fllFILT:! 

IFIOBSI  TINCONl.LT.. 001 IIFILTSI 


c 


IUFUTSR 

r  IOOSI  TJNCON.S.I.LT  ..1 IIUF  ^T  sf 


c 


ificptr.lt. .oiicptrsi 


wri  if  i*  ,*oso  innser,nplins  ,cptr 

»oso  fimMiM'  mm  uwiiict  - - 

«  *i  no.  or  riEiuiK  limits  s«,i*// 

•  •  CONDUCTIVITY  ;tu  fACTCO  s  *,612.k//> 

c 

(  M«*n  PRESSURE /CYCLE  limii  carosss 

e 

c  rnr,mr,Rii#iiiiia - 

c 

C  MOM  /MOT  m  UPPER/LOWIR  N[SIUl[  LIMITS  TO  OE 

C  APPLIED  TO  INPUT  OATA  CYCLES  OCTnCEM 

C  CYCLES  Ii  m  10.  EACH  SEQUENCE  OF 

C  CYCLES  SPANIINO  THIS  PRESSURE 

C  IMTEMvAL  IS  MAOC  INTO  AM  OUTPUT  OATA 

C  SET  00  POqFIlE.  (PRESSURE  LIMITS 

C  |  YE  IRRELEVANT  POP  RCMOPsi.l 

C  SELIAP  IS  A  SLOP  FACTOO  (USUALLY  AOOUT  2  Oil 

C  WHICH  IS  EIPLOVCO  TO  POEVCRT  SMALL  PRESSURE  ptUC- 

c  tuatimmS  Pnom  toippimo  up  the  pocssuoe  limit 

C  LOOiC.  I  T  INSURES  THAT  A  DOWNCAST  HAS  EhOLO 

C  SEP  ONE  THE  SUOSEQUCNT  UPCAST  IS  OCCOsNjECOi  ETC. 

C 

c  11.12  THEME  IS  A  OLOOAL  CYCLE  I  NOES,  SAY  lE.  WHICH 

C  IS  SET  TO  I  AT  THE  OESIMMlNS  OP  AN  INPUT  OATA 

C  SET.  IT  IS  INCOEMCNTEO  POO  EACH  CYCLE  CnCouNTCOEO 

C  WITHIN  THE  DATA  SET.  (AM  INPUT  DATA  SET  TYPICALLY 

C  SPANS  NO  IE  I  HAN  ONE  SEOHENT  AND  OCCASIONALLY 

C  HONE  Than  CNE  pile. I  THE  fiost  POOPIlE 

C  START,  OEFWEo  NT  PT0P,  |S  NOS  SOUCMT 

C  UNTIL  IE  IS  SE  II.  AFTCO  EACH  POOPILE 

c  is  completed,  tp  ie  is  se  u,  the  meat  poessuoe 

C  LIMIT  CAN!  JS  INVOKED.  WHEN  THE  CA*OS  ArE  E»- 

C  HAUSTEO  00  INPUT  DATA  SET  ENDS,  THE  MEAT 

C  OATA  SET  CARO  IS  READ.  IP  THE  ENTIRE  INPUT 

C  DATA  SET  IS  TO  RE  SCANNED,  SET  !UI2=0. 


11:1 

CO  JO  L: I.NPLIMS 

RE  ADI  S,  10001  (PL  1 1,11  l,Xsl,Sl,(  IPLIRII.II  I ,  Is  1,2 1 

IF  1 1PL1MI I  ,1 1 1  .LI  .01 1PL1M  1 1  ,11  >S1 

IF  |  IPL  1M(  2 ,  il  I  .IE  .01 IPLIN  (2 ,11 CIOOOOOO 

kRIlE|N.A010INHSER,|PLII.IIl,ISI,SI,HPLIN|l,III,I=1.2| 

•0 10  FORMAT! |K, ’SERIES •,( 1,12, •( M I, *PTOP=.,PS.O, I *, *1*. 11, *PROT  =  *, 

•  FS.O  ,*  (  *,!*, *DELT»P-'> 

CF$.0.*( *, 1A, 'CYCLE  1? * , 1 T, 1 S,*CvCLC2:*« IT  I 
C 

11:11*1 
TC  C0NI1NUC 
C 

irilFlLT.EO.R.ANO.lSW.EO.ll  *0  TO  VS 
SO  TO  (2|, 22,21,1001, 1PILT 
C 

C  MHO 1  FILTER  SENERATISN  SECTION 

C 

52  1FIaWTT.LT.1INWTT=1 

CALL  SE«WMFINNTT,T!HC0N,NT> 

C 

WRITE IA.A0T5INWTT ,TIMCON 

oois  formats/*  wooes  hole  deconvolution  chosen*// 

•  *  NO.  OF  NEIOHTS  =  *,!*, 

•  «  TIME  CONSTANT  S  *,6n.S//l 

c 

NT:(NWTT<1 1/2 
PC:R1 
AW!C:nwtT 
00  20  i:i,NHTC 
50  WCIIU1./NWTC 
C 

SO  TO  21 

e 

C  SPECIAL  filter  input  section 

c 

21  WRITE  |0,2110IIUPlLT 

21(0  FCMNATI//*  SUPPLIED  TEMPERA  MR E.  CONOUCRIVITY  f  IlTcRS* 

•  /*  TO  RE  READ  ON  IN  IT*  .1  W//  I 
C 

C  ••••••  SPECIAL  riLTERSl  NO.  TEMP  UTS,  NO.  COHORTS, 


tw  ui«  tcne  ttti  cmo  mt. 

THE  SC  VALUES  AM  ICAO  EMU  CAMS  IE  Tl«M:4. 
OTNEAHItC  II. C.  .1 

THEY  A  AC  A  (AO  EAM  IM|T  A.  IE  TIMM  AE  -I. 
INIS  AEAA  STATEN  CUT  IS  AM  IE  ECO. 

ACAOI  lUEtt  T.lOQO.CNOtf  VVI  NUTT  ,NMTC  .LAST,  INI  |I  1. Ixl  ,NNTT»  , 

•  UCI  It.Pl.NMTCI 

Mill  (A .7401  INN TT tNH TC.LaAT 
TACt  fOANAfl/*  NO.  OE  TENECAATU  AC  VCIANTSl  •»!  *// 

•  •  NO.  OE  CONDUCTIVITY  KS*1S>  Mj// 

A  •  TCNECAATUAC  LAO  :  ',11,'  0*1*  CYCLES*/! 


XTMNMTT.l  1/2 
NCtINH |C>I 1/2 
«ICM»:NT*  2-| 

ACCHE:nC*2-1 

IEInUTT  *£q .NICHE .ANO.NHTC .CQ.NCCNNI  SO  TO  n 

c 

AAIIE|A. 20001 

21(0  EORNtTI//*  NO.  Of  UCIANTS  NUST  AE  000*  TOY  *A*IN.*//I 
S10E 
C 

C  C«£* It  EiESSUAC  SNOOTNINS  EILIEB. 

C 

<1  Ail  IE  I A  .2 100  ) 

21  CO  foanAll/*  ENESSUAE  SNOOTNINS  EILTCA*/| 

C 

CALL  AE* LIINE.JA.JA.JC.E.MEI 

c 

«U1E:2«NE-| 

SURE :0. 

Sue  c  to. 

SLNT-O. 

Do  2!  1  =  1. NOTE 
as  sun e t June. ue i s s 

DO  2T  It  I, NUTT 
21  SUN1:SUNI.UT|II 

00  as  Itl.NUTC 
as  suNCtsuNCucm 
c 

Ail  IE  IA.2TQ0INMTEiSUNE.NUTT  iSUNT.NnTC.SUNC 
2TC0  (OSNATI//*  SUN  OE  *,IS«*  ENCSSUAC  HE  1AN|  J  =  Mll.M/ 

*  •  SUN  OE  •,!!,«  TENECAaTUNC  HEDnIS  X  >,011.0// 

•  •  SUN  Of  CONQUCTIVITY  hCIANTS  S  ..SII.A//I 

c 

c  #•••••  EAOCCS  ONE  OATS  SET.  AAAAAA 

( 

ICO  CALL  ACCC sSlN.LNl 
ISHX2 


(0  10  ts 


YYV  Hi  I  IE  lA.AOOOl 

AO 00  f OiNAT 111. * AAAENO  07  JOOAAA*! 

STOE 

C 

YS  Ml  IE  IA.VSOO! 

VSCO  fOiNATI//  •  EiEYIOUS  7ILTC  A  TECCtflEO  ON  INITIAL  EASS.* 
•  •  INCONSISTENT. •//! 

SIDE 


AEiT.S  EfAUOTAEESO. ACCESS 


PrtUOTtPNOOl 

1 

2 
s 

4 

s 

1 

* 

10 

It 

12 

II 

1« 

IS 

IS 

IT 

IS 

l« 

20 

21 

22 

21 

24 

2* 

2* 

IT 

2* 

2f 

10 

11 

12 
11 

14 

15 
U 
IT 
It 
1* 

40 

41 

42 
41 

44 

45 
U 
4T 
40 
4* 

50 

51 

52 
SI 

54 

55 
St 
ST 
St 
St 


II. ACCESS  _ _  _  . 

COMP  ItCNIB  110*21 
SUOMOUTINE  UttllltfUl 

C  COMMON /UNA  1 /NINO, IMO  1 201. AC  MOP  .  NPLINS, PCI  2.100 1 « IPC1NI2  >1001 , 

tHMM*l.tmiTl20l>ICtLU 
CONNON/ROOCP/FOOCRItOt  _ 

cohnon/rooci/ioocriooi 

COMMON /0§0  CO/ AQOCNI 1001 

COMMON  /  MMOO  /  LM.NM.MRM  «NWM  ,NNFM.NpM.N|N.NAM,lPMI20l 
COMMON  /  ONOO  /  tO,NO,NOO,MWO,N;?rN 

COMMON  /MOOCr/TOoCUltOI  /MOOCI /TDOCMI 401  /*»CA/A0OCU|10O» 
COMMON  /  MOtTt  /  VMttOOIS  I  _  _ 


COMMON  /  0100S  /  NSSO.NSOM.WNO.NNNMtlMlO.NNr.NNI  JMA, 
DOTS  MMNO.NNNM.NNIO.NNr .NNI MN 0/1 000. 1000. 12,40.40, T42 
GOTO  NrM.NlM.NAM/40.40, 100/ 


INST, IMS! 


C1MCNSI0N  M 02001. ICT 12001  ,RCc0UT|4,200I.RCMUT  1201 
ClMCkSlON  01201, S0M0CCI4I 
OOlO  JO  ,J0  ,JC/2.1  ,1/ 

LOB  1C0L  IlOLU 
MCMOr ISMCNOO 

100Cb<SI=0  _ 

iriN^D-l  lLN=l  _  ' 

IMKCM0P..E.1IS0  TO  >00 
k=0 
1 10  -2 

IOC  CONIINUC 
1D0CJ  1=0 
t»TZ=0 
H0CMI4IS0 

Sr 0 NMSM .HE .NNIhO IIOOCM I  SI =0 
100CMI5I  =  100CMI  Sl*l 
jP*=2 
4020*0 

IT  tINOltl.EO.OtSO  TO  J2S 
CO  124  I|=l,t 

124  1M1mOUI>.NC.OIJn20=JN20*1 
)2S  CONTINUE 
NM=0 
IM=NIN0 

IMMMOi.SI.NNNMINMMOI=NNNM 

00  lot  J=1,LM 

1PmUI=IPNSCTIjI 
SOI  CONTINUE 
lNII-1 

ISO?  _ 

If  IN.LI.1I  JSM=1 
1E0LIC  .FALSE. 

1CI  to  10  H.2I.ISN 

C  SlR01SHl*lHROUSM  PNOCCSSINS  SECTION 

1  COLL  CTRCClStO.O.lUOl 
RUO:lUOo2*l 

C  CEN lit  OPTION 

TO  irilEOLUlNETUNN 

COLL  OCRl  VE  IQ.lNl  T I 
lkll=0 
(0  kk=NN«l 

IJ0=INm-1I4LM 
CO  II  J=l.lM 
14=1J0*J 
II  ON  II JI=Ql J I 

IFINN.EO.NMNAO.il  00  To  PS 
tO  TO  I 

C  EEC  IM  HON  SECTION 
C  SET  riRST  RECONO 

C  lbO  =  l  t  UP  MOTILE  IUOsjO  00 MN  PROPILE 

2  CALL  RETRECItVQ.0. 1V01 
OUO  =  1UO*2-1 

C  !Ct  FIRST  INTERVAL 
P=Q I JPOI 
IP:P*IO.*.S 
1PUIP 

IPO  :IPI*LN*N4RU0 
IP2  =  IP0*LN«40U0 
IPP=1P 

out 

MrlNON 

RMT=2*N*I 

IFlfc.EO.il  St  TO  T2 

IF  IN.EO.NPOEV.ONO.LN  lEt.lIPNT  I  SO  TO  T2 


159 


I) 

H 

•  S 
It 
IT 
■I 

•  I 

50 

51 
12 
IS 
II 
IS 
II 
IT 
II 
II 

100 
101 
102 
101 
101 
10S 
101 
10T 
101 
101 
110 
111 
112 
US 
111 
US 
lit 
UT 
111 
111 
120 
IZ| 
122 
121 
121 
12$ 
12t 
121 
121 
121 
ISO 
111 
1  S2 
US 
111 
1SS 
lit 
1ST 
IS! 

111 
110 
111 

112 
i*: 
in 
us 
us 

UT 

HI 

HI 

150 

151 

152 
1SS 
ISI 
ISS 
ISI 
1ST 
IS* 
IS* 

ito 

111 

112 

its 

111 

IIS 


FOU.Sn 
SI  UCII.ISIII 

isii  fiutii/'  occiMTioM  ruuivi 

CALI  tlNFLllH,J|,J|,  JC.rO.VI 

limit  .non 
tici  roia»Ti//» 

IPR£V:N 
LtPICV:LN 
CO  1C  12 

C  UT  It  ll  RCCORO 

TS  CALL  S£IRCCItlO,Q,IUO| 

100:100  *2-1 
IS  P=Q»JP«1 
lp=P*IO.«.S 

C  CMCCK  FOR  P-REVCRSAL 

T2  IMIlP-lPPHKuO. Ll.Ol  CO  TO  TS 
IPPSIP 

C  CMCCK  FOR  P  OUT  OF  1NTCRVAL 

1FII1P-1P1IIKU0.LE.0I  CO  TO  TS 
1FIUP-IP2HKU0.SE.01  SO  10  CO 
C  IOA0  RCCORO  INTO  FILTER  ARRAT 
UP:  I  IP-  IP  I IIKUO 
IC]|JPI:ICI  |JP |*i 
CO  It  J:S,I 

Tt  SEC  OUT  (  J,  UP  irRECOUT I J.JPI *Q  Ul 
RCCOUTU.  JPURECOUTI  l.JPI  >0111 
CO  TO  TS 

c  complete  averaging 
(0  CO  tS  K=K1,  NmT 
KCl:lCTtHl 

lrlKCT.NE.Ol  SO  TO  It! 

KKIH0L:KNTM0L«1 
KN12=KNT2*i 
00  Its  J:s.l 

ICS  GECOUT I JtK |:SAVRCC1J I 
RECCUTt l.K I-SAVRECI 1 1 
(0  TO  65 
1C1  CONTINUE 

IF  1KNT2.CT  .1  lURITC  It  ,SII1  IK  NT  2  ,P0 
Sill  FCRHATtll.'GAP  0F*,IS,1(,*AT  P:*,C12.t| 
KNT2-0 
CO  IS  J=S . I 

SCCOUTI  J,K  |  iRCCOUT  t  J  ,K  1/KCT 
t!  SA  vrECIJURECOU  Tl  J.K  I 

RLCOUTC  1,  Hl^RCCOUTll.KI/KCT 
SAVRECI1URCC0UTII.NI 
AS  C0NI1NUC 
C  IF  NO  FILTER,  SKIP 

IFIN.CQ.il  GO  TO  tl 
C  CON  To,  VE  FlLlEN 
KSUM:0. 

«C  1:0 

OC  tl  j:s,l 

tl  KtKcUTIJUG. 

RLKOUT 1 1 1  —  0 • 

CO  tt  K:i,NVT 

IF  1 ICTIK | .EO.0 1  SO  TO  tt 

vsun:h  sum  *  U  IK  I 

.Cl:JCT*i 

CO  IT  J: s, I 

fT  REKOUTIJI:RCKOUT | JI*RCCOUT| J.K IINIKI 
RCROUT II I =rEKOUT 1 1 1 *RCCOU  T 1 l.K l*H IK  1 
tt  CCNIINUC 

IF  I  JCT.CO.O.OR.JCT.CQ.NIIT  I  GO  TO  tl 
CO  tl  J-S • I 

tl  REKoUT1.1:rCKOUT*JUNSUN 
R(K0UTUI:REk0UTIiI/HSUN 
EC  10  tl 
tl  CO  t2  £S,I 
t2  RE  ROUT 1J|:RCC0UT I J,| I 
IEK0UiUI=RCC0UTII,ll 
tl  CONTINUE 

FC:FLOaU1PO»/10. 

RCKOUT  I2I:P0 
C  RCCORO  NON  complete 
C  CIR1VI  OP  I  |0N 

CALL  0CR1VE  IRCKOUT  ,IN1T  I 
IN  I  HO 


11  Kv:nn«I 

1 JC- 1NV-| IILU 
00  SI  J:|,U 
I  j:  1  JO*  J 


C 


it* 
iti 
IM 
U« 
1  10 
111 
Hi 
111 
IM 

ns 
in 
I  IT 
lit 

uo 

III 

III 

l«J 

IM 

ns 

IM 

III 

IM 

M* 

ISO 

111 

III 

111 

III 

ns 

m 

in 

ns 

i« 

200 

201 

202 

201 

20* 

20S 

2c* 

20T 

20* 

20* 

210 

211 

212 

211 

21* 

215 

21* 

*1* 

21* 

*1* 

*20 

221 

222 

221 

22* 

22S 

22* 

221 

22* 

22* 

210 

*J| 

212 

211 

21* 

21S 

21* 

211 

21* 

21* 


2*2 

2*1 

2*« 

2*1 

2«* 

2*1 

2«* 


si  iiiunHiwiui 

1FINN.C0.NN**I«I|M  II  *s 
C  •($((  »tl  SUMS  *M  MIMICIS 
•  (  |*fl=  l*C*M*»UO 
tM|:lM0-CM««M|lO 
I*2=l*0«lR*N*Ru0 
■  illllKI 
*Nl*M:NNt-N 
0*  «S  *=l  .NN1NN 
**:*«* 

1(1  I*|:IC1  |UI 

■ccauf  it  • «  MKcurii ,« » 

(C  «*  j:|* 

«S  *C£0U1|J,«  IMCCOUtlJ.SNl 

rc  «»  ci  i« 

•*:*Ml>**l 
let  is* I=q 
HC«l||i«ltUI. 

00  «*  J:l,* 

«*  *(COUttJ(**t=Q. 

*1  IMUO-ul  KIM* 

iruoocui.es. use  i*  too 

(0  10  IS 

<s  cm  i(ncciiioiiiiu«i 

MUQ-ftf  0*2*1 
**  ioocni it=o 
c*«  no *i  siurr  hc*c  •• 

«T  CON  lINUf 

c  sMiri  n**2 

iruN2k.eo.oiso  to  122 

IM  I -HU- 1 

jmfi=ju2* 

00  |2I  J=l.NNHl 

«*n*«r  ti=MMj*«ri«LHi 

.BVF1=J*«F l«LU 

121  C0N11NUC 

122  CONIINUC 

ir  INN.  SI  .NNNSi  |NN  =NUH*S 
NN*  N-NHB* 

*HFN=NhFR 
CO  10  1*1  .*21  tlUO 
*2  100Cn|211=0 

IOOCnI* I-IOOCu |« I* 1 
DO  9*1  10 -l  »*0 
*t  1  foocn«Mi:Foocriioi 

1CCCNI2I=IG0C*I2» 

ioocniii:iooC*i 11 

Io0c*ll0l=l00c*liol 

loocNiiiisioocmiii 

lOoCN I 12I-100CR 1121 
CO  *2  10-1, N«* 

«*;  kOoCNiiqi=MlOcRliai 
CkLL  2jRI1|2.1F,0| 

NNRMP-NNO* 

so  10  1*S 
51  100CH|2I|=I 

1DOChI*I:IOOChi*I«i 
co  «*]  10  =  1, *0 
9*1  FOQCU  I  lQl  rFOOCR  I  IQ  I 
10  (CM  1 2 10 1 00  CRl  21 
ioqcui  mtooCRi  si 
loocMitomoocRiioi 
100CN 1 11  IstOOCRt  1 1 1 
IOOCNI 121 OIOOCRI 121 

co  ««  iqoi'NM 
*•«  tOOcUlIQI  :«OOCR«|OI 
C*LL  2j  R1 1  (ItlFfOl 
*NBN*=NH*U 

1*!  ir  INNINOL.sI. OIHR licit »**10  MMINOl.nSH 

**10  FQRN*11U.ISiU.*M0LCS  IN  0«l*  RCPLtCCD  MTh  MCCCOINS* 
(  1K«*r€C0R0  IN  SCSNCRI  *tII,ISI 
RNInOLSO 

c 

««  »N=l 

IOOCNI  1 1=0 

c  shift  *cco«o  nu*i  10  1 
»NR**L=MIN«**tM 

co  ia  j=i.iu 

«NIJ1:«M|NNN*1L*J1 

12c  cmiinuC 

in  IF  .61. 1 1  so  10  *0 
SO  to  llO.tll.lSU 

to  coaiiNuc 


2N0 

100CMI  1  IS! 

250 

IOOC*t=l 

SSI 

1*11=1 

2S2 

(0  *00  lH*«s|tNMT 

251 

IC1  IIK*tt=C 

2SS 

to  OS  IRSsI.S 

2SS 

SOS 

SECOUT 4  INN  ,IM*9  ISO 

2SS 

IF  4nM.Co.oISO  TO  NT 

2ST 

SO  10  tT 

2*0 

cs 

SS11  {IS tOOOOl  IF 

25* 

toco 

FORM*!  1// *  PROBLEM  IN  2MRITV/I 

2*0 

STOP 

2*1 

C»l 

sprt.s  rrMonr(H.((i>cc 


rMuoi«num<uiKC 

1  C  CMPHERI0I»S  =  1I 

2  SUBROUTINE  SET»CCC*«0.lUO» 

i  c 

*  DIMENSION  01201 

S  COMMON/ tfH*  f/NINOtSNO  1 201  t«CMOP  (NPLXNS«PLI 1,100 1 t  IPIINC  2  ,1001 

*  •  NNHSX.  IPWSET 120  1 1 ICOLU.N  MER 

T  COMMON /CtL  •  /*  .NT.MC.lSS  It  MM  *01  .MTISQI  |UCI  *0 1  .PDIF, 

*  •  IDIF  .COIF.cFTR 

«  CGNNON/RHOR/LR.NR.NOR.NNBR.NIFR.NFR  ,NI*,N*R(1PRI20I 

10  COMMON  /MOST*/  VRl*. 11001 

11  COMRON/IOOCF  /roOCOIMQI 

|2  CONHON/ROOCX/lOOCNCMOl 

|1  COMMON/MOOC«/«OOCRIIOOI 

IN  COMMON  /Ol*  SS/mSGR  tHSSM  ,NNNN  .NNNM.NNIP.NNF  ,N»I,  NN«  ,IR$T  ■  IMS1 

IS  CM*  NNNR.NNIP.NNF.NNI.NNt/BOO.IO.NO.NO.lOO/ 

U  C»l«  Uu/l/.JM.JP.JT.JC/l.za,*/ 

IT  0*1*  1IU/*/ 

10  CM*  ISEtf  INCXI/Itl/ 

I*  easiest,  l eolu 

20  CO  TO  l*,7|,ISCT 

21  *  CCNlINuC 

22  lOOcMI 1  ISO 

21  11I:MT«L*SI 

2*  ll  =  MC 

25  IF41lT.ST.MCIX|=IIT 

2*  1F4HP.GT. MIIIPIP 

2T  MMU2NMI-I 

20  MHC:2*NC-I 

20  MHM -2NMP- 1 

10  ISC1=2 

11  I»SI=2tI|-| 

12  T  CONIINUC 

11  C 

J*  C  1*11 1SII2C 

IS  C 

IS  CO  10  M,2|tlSM 

IT  I  *Sm:N 

10  CALL  IUSET 

It  TOC  SMSERUNIl  |«1 

no  ifinmseoi.co.nmscrioo  to  toi 

N I  CALL  2*C*0IIIU.IF,0) 

N2  IF  I  IF.  CO.  I  ISO  To  702 

Nl  IF|100cMIII.C0.IINl0l:N10I*I 

NN  CO  TO  TOO 

NS  TO!  CALL  111  SET 

NS  CO  TO  TO 

NT  TO  I  CONIINUC 

NO  JDOCRIIKO 

NO  ISM-2 

50  *=I»-1 

51  1CPF:Q 

52  I0L1=-I 

SI  SO  1C  S 

SN  N  MSMSl 

OS  TO  1  IIsIC-l 


SMiV* 


c 


St 
ST 
St 
4* 
tO 
tl 
t2 
tl 
t« 
tS 
tt 
t  T 
tt 
t« 
TO 
Tl 
12 
n 
Tt 
1$ 
It 
TT 
Tt 
T  4 

to 

tl 

IT 

tS 

tt 

tS 

tt 

ST 

tt 

S9 

to 

tl 

t2 

tl 

tt 

tS 

tt 

tl 

tt 

tt 

100 

10J 

102 

101 

lOt 

IOS 

lot 

10  T 
lOt 
lOt 
110 
111 
112 

11  3 
lit 
IIS 
lit 
111 
lit 
lit 
120 
121 
122 
121 
12* 
124 
12* 
12T 
1  2  » 
I2t 
110 
111 
112 
Ill 
lit 
US 
lit 
111 
lit 


s 


CtlL  lUSCT 
COM  I  IN  UC 
MNNR1=MNNR 


2  11- It  ♦  I 

IiL:HIC*ll-2*»CPr  1/10001*1 

IFIltl.CQ.  lltl  ISO  TO  ItO 
IMI00CRIII.Cfl.ll60  TO  20 
IMIoLl.cO.-llSO  TO  1J1 
1041:2*11-1 
irill.Cfl.H60  TO  111 
1MC6=NR*1 
to  $  I :N0*I6 ST-1 
-  -  0 

co  in  mtcs.MtsT 

j:j*1 

CO  1 12  11=1. CO 
I  32  «0I1I,JI:  t0ill.ll 
I  11  Cut  UMUC 
131  C0MI1NUC 

ull  Toctotiiu.ir.oi 
irur.cci iso  to  rti 
»rsi=nr«ios  i-i 

JlM= IN (XT  -  1 

If  I  j)*.tO.OIJlt-l 

If  IIC.LT.  IPLIMI.JINIIOO  TO  ItO 

1RS1P=IR4I*IPlIhI1.J1|  1*1 

.  J  :  u 

CO  Itl  J: IOSI, MR  ST 

-<•  1 1 
-PI  :J*I 

J*i  N  i:  J —  | 

ir non . co.-i. too. JNi.tc.iosir iso  to  iti 

0irp'l6sl«Rl2.J>-*OI 2,JMI II 
ci r latest vRii.Ji-vRii.JMi H 
circ:tasi«Rit.Ji-«ott.jNi t» 

ClFPP-ttSl  MR  12,  JP1  |-tO  t2.  JN1I1 
C1FIT:*BS1«R11,JP1  l-tR«l.  JN1»1 

OlFLC-tBSItRlt.jPil-tOlt.JNllI 

iriciFp.LT.POir.oR.POir.LC.oi60  to  it2 
lFIOirpp.se. IP0IFII60  10  1 12 


tl 


MR1  ICIt.tlt2lM0R.JU.  IP0l2l.ltlt2,Jl.ll0l  2 ,JH1I 

i2  roflNrnu.'set'.it.is.i*, •cKLC'.is.u.tt.w.r 


l»  2 


It 


Itl 

ltc 


l  'REPitCCO  tT  *»l*.r  10.1,11 12, *s  *1 
to  1  i  i  JUtO  1 2. UNI  I 

irioiu.LT.Toir.iR.TOir.LC. quo  to  its 

IF  iQlFlT.iC.lTDIFDEO  TO  |tl 

«UUIt,tIt2|MtO, JJ.IPOlll.TOI  J.JI.VOII.JNII 

tRl  l.JI-yO  1  l.JItl  I 

ir tcirc. d -coir. or. coif. lc.oiso  to  jti 
if  ioircc.ee. icoir  H60  to  itl 

MR  lie  It,  tltTINOR,  JJ.lPRlt  |,«l  t  «<!!,**  It  .Jit  II 
IR  tt ,J  UtR  It,  JH1I 
CONllNUC 
CONIIMUC 

ifuoocrii  ».co.iiNioi=Nioi*i 

iumr.lt.  ioooiicpr=iooo  jio.icpf 

lien 

ir  i  itti.co.-iMCcui  *i»$t-i 

MMNRi:NR 
181  1-1  (L 

iuir.co.oi60  to  ito 

IUir.6T.ll60  TO  22 
CO  10  t 
CONllNUC 

MRticit.toooiir 
toco  FORpSTtll,  ••••CRROR  IN  2RCtO,  lF=»,II,*t*t  I 
SIOP 


10. 1.  It. 


22 


ItC  CONllNUC 

c 

c  c cnvolmc  filtcos  mitm  otit 
c 

CO  II  j:l.t 
II  0 1 J  1:0. 

ell  I -MRIJM , ICC  1 
C0  12  1=1 ,NNT 
1T:ICC-NT*1*L»6T 
12  ClI|:0|ll«M«lJT,lTltMTIII 
1C  II  1=1. NNC 
IC  =  1CC-HC*1 


22 


1 J* 
1*0 
141 
1*2 
1*1 
1** 
1*5 
I** 
1*7 
1*0 
1*4 

no 

151 

152 
151 
IS* 
155 
15* 
157 
15« 
154 
1*0 
1*1 
162 
1*1 
1*4 
1*5 
1*« 
1*7 
l*i 
1*4 

170 

171 

172 
171 
17* 
115 
17* 
177 
171 
174 
110 
1*1 
1*2 
1*1 
IS* 
1*5 
It* 
1*7 
l*t 
1*4 
1*0 

141 

142 
141 
14* 
145 
14* 
147 
14* 
144 
200 
201 
202 
502 
20* 
205 
20* 
207 
20* 
204 
210 
211 
512 
211 

214 

215 
21* 
217 
21* 
214 
220 
221 


Cl*  l=gl*l«x*IJC«lCI*NCIII 
11  commit 

co  i*  1  -i  ,  DM* 

IP  =  I£t«l-RP 
«i2i=o» 

1*  CCKIINUE 

c 

c  *mi  mnufacturers  cell  connection  and  crin 
c 

CALL  N»CCCIQ42|,a«Jl,04«l,CFIRt 
C 

ICE-ICC  •! 

C 

C  *  CMO 111  I  ALL  OATA,  KCHOPrljUP  *N0  COHN  CASI 
C  N  CHO 1  - 1  DOWN  CASI  ON  Li  KCHOP  =*  lllf  CAST  ONLY 
EO  10  lit, 100, 100, 1001, KCHOP 
C 

II  IFIlNExT.GI.NPLlHSlSO  10  200 
It  1  -IPL1H1  i.InEYH 
1ENEYT-1PLIN12,  INC  AT  I 

lf«  IE. LI, IE  1160  10  2 
1M1E.L1.IENEI  11*0  TO  25 
IN£«7  =  In£*T«1 
*0  10  14 

c  e»o  or  lu  on  file 
c 

20 C  I*£xl=l 
20  ISri-l 
1SEIU 

lc  OLOr.TRUE. 

A101P1 :N10 1*1 

irilDOCRliI.EQ.llNIOlPUNIol 

lFiNHSCR.NC.NlDlFllMRI1Elt,4*S2tNIDlPl,NHSE* 

4*12  lORHAIl///  IX ,*CONGRA  TUL*T!0p6  TURKEY !  !  I  TOj  HA  VC  TtlEO  70  PROCESS 
t  stnlES*, 11, 'INSTEAD  OF  SER  ES  «  ,1 1, 1*/'  CHECK  PRESSURE  LIHiTS'///l 
tttliRNl 
25  RE  TURN 
C 

ICO  FrOIJPl 

c 

to  10  4110, 110, U0,122I,KSH 

c 

C  SET  UP  P  -LI Hi  IS 

c 

132  FltiPLll.ll 

IF  4PIT.lt. C.lPlTrAINUPM  .1 

P2t=PlI'PLIl,ll 

Flt:PL(2, 11 

IF  «P1R.U.0.IP1R=*4000. 

P2R:P1R-PL<1,1I 
ItlrlPUHU.il 
IE  REX  t=IPLlH|2,ll 
IN  Ex  I  -2 

IF  lit,  Lt.IElltO  TO  2 
C 

C  SEARCH  FOR  FIRST  PROFILE  START 
11*  CO  TO  11*0,1*0.  1*0.1*11. KCHOP 
1*0  IF  IP.lt.P)  1|60  TO  117 

to  TO  U*l. It  1,1*2, 1*11. KCHOP 
l«|  IF  IP.tl.PlRIGO  To  11* 

1*2  CO  10  2 
C 

C  SKIP  FOLLUING  RUBtlSH  IF  PL1N  EWAUSTtO 

c 

c 

c 

nr  ifiie.lt. ieiioo  10  2 
c 

to  10  111  I  ,11  7,11*,  I  1*1, KSK 
c 
c 

c  CHECK  FOR  P  INSIDE  LIHITS 
C 

III  lFtP.6T.PIBt  60  TO  US 

iup.it. phi  60  to  in 
CC  To  U!*,USi,IuO 
t 

C  t»0  OF  PROFILE 
C 

112  *S»:2 

to  to  1151  ,1S|,2,1S1 I, KCHOP 
I  IS  KSK  rj 

60  10  11*1, 1S1, 151,2 1 , KcHOP 
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222 

1S1 

COkllNuC 

223 

IC01U  =  . FALSE. 

22* 

if  ur.  st.inirtz.MrLiNsiiM  to 

225 

iMIOOCRTtl.EQ.IICo  70  20 

225 

irilC.tC.lCNCRIIsO  TO  124 

227 

iruNtn.sT^ritMSiso  to  124 

22t 

c 

iscrenent  rims 

229 

fii=ri» l.Tkfiii 

230 

F)S  =  F>1 12,INE1T| 

231 

p2 i=pi  t*pl  1 1. 1  nc  st  i 

232 

F2(:P|6-PL t J.INClTl 

233 

ill  :inimi.i«if  ti 

23* 

ICNest  =  iPi  In! 2  •  Inert  • 

235 

INEAlzINC  RT*1 

234 

124 

COHTINUC 

237 

RETURN  1 

239 

c 

237 

c 

si 

TCP 

2*0 

c 

2*1 

117 

NSu  =  2 

2*2 

If  IP  .1  t*P2l  1  CO  TO  2 

2*3 

1K0-2 

2** 

lit:) 

2*5 

1  IS 

CO  10  |2S«2St2S,2l,RcM0P 

2*4 

c 

2*  i 

c 

St 

SOI ION 

2*4 

c 

2*9 

1  is 

NSN  : 3 

250 

If  (P.SI.P20I  60  TO  2 

251 

luo-i 

252 

NSN  =  1 

253 

1  u 

to  10  <2S,2S,2,2SI,NCM0P 

25* 

c 

255 

EnO 

•PRT.S  Pf*U074pS6D. lusty 


PF*U074PR6D<ll.lUScI 


1 

CONPllER|OISC:il 

2 

SURROUllNC  IUSET 

1 

CON  PON  /10f  HE  jNlNf  TLtlfEfli  }20  1 

4 

OInCnSICM  ISRltSI 

S 

csis  imcfo/ 

t 

If  Ut  =  lfILC«l 

7 

I  s  « If  IlE-l>*6«l 

S 

111-1*3 

9 

If  tlfllC.CI.NINFHICO  TO  900 

10 

ENCODE Sl6t  12.LSRl»IIFER|J  1,0=1  til  11 

11 

12 

FORNST  1  *SU  SC  4  •  t  *4  S6  t  *  .  *1 

12 

Call  ERTRs NIS, * SfREC  a S  4.  .  •  1 

11 

CSll  CRTRSNtSalSRII 

14 

kKltt<A«tOOOI<IflOIJI,J=I ,1111 

IS 

60C0 

fCRRAT 111, •FILE  SSSI6NE0*  ,1I,4S4I 

IS 

SC  URN 

17 

9  to 

If <1FIIE.C0.11  60  TO  901 

IS 

1  =  1-4 

14 

111  =  111-4 

20 

SRIlEIS.ROCOllIfCRIK 1,6=1 ,1111 

21 

90  CO 

|ORN  ST  <//  •  PR06RSN  TERNINS1ED  OFF  ORE  RESCHIN6 

22 

»  •  ST  CNO  Or  INPUT  PILE  ',4**1 

21 

SIOP 

24 

PCI 

,KI IE  IS ,9001 1 

2$ 

90C1 

f  ORNSt  <//  •  NO  INPUT  rilE  NSNCS  SUPPLIED.*// 

26 

f 

•  '  Itpul  Flic  ASSURED  ON  IN I  T  4*//t 

27 

Rt  turn 

28 

END 

<m,s  pmuotsprso.sennhf 
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« 


I 

I 

I 

J 


SUBROUTINE  •CNHHMHMT.IIHCOM.M  I 

c 

C«*  -  WOODS  MOLE  FILTER  S  lit  , 

C««  IIMCCN  -  I INC  CONST  OUT/ ft  NPLlNO  TIME  INTERNAL  00  INSTRUMENT  . 
C**  W  :  WOODS  HOLE  riLIER  U*tl  NS  HE 

c 

DIMENSION  Wilt 
fAwt=FLC4T<NMT| 
lfl.Wt.CQ. Cl  OETUNN 
CO  10  |K  ,NWT 

r:rLO«mi 

t  ii  ini*  FNwT>«iiTiNC0Mi*ii2.*r-«.*(rNwT*i.i  1/ 

I  IFWWlWlfSWJWW2.-l.lll 
10  CONIINUC 
REIUON 

two 


RFRT.S  OfwilOJWp  f  40  .OENf L  T 


l 


1 


! 


1 


Pf  wUOlSPRGD  t IT «  f (Hf L I 


CWWWWW  WWW  WWW#  WWW  WWWWWWWW**  *WWWW  WWWWWWWWWWWWWWWWWW*  WWW***** 

c 

C  SINPLI  IS  C  SSE  NTI  ALLY  SENE  Ml  OHM  |N[  rE»»SA  PACsSOE. 

C  I  SEE  1  f  AST  AND  EAST  11nC  SCOIES  ANALYSIS  SI  NCSU  - 
C  E  CITED  BY  0.  A.  RROOHS  -  AN  NCSU  INTERNSL  DEPORT  I 
C 

»  -  WALf-WIO  TM  Of  OILIER 


.A  -  l.-l  fOR  onE-SIDEO  FILTER  fSOTNS  TO  THE  RIGHT, LEFT 

2  ros  a  symmetrical  filter  or  izm-ii  terms. 
wi  •  0  fOR  las  WINDOW  application,  lsrsesi  TERM!!. 

1  foR  LOW  PASS  FILTER. 

2  FOR  MIS H  PASS  filter. 

.c  •  1  rartlett  taper 

2  PAR2EW  TAPER 

3  COSInE  itureyi  taper 
N  LA  NC20S  TAPER 

S  LAWC20S  SQUARED  TAPER 

flLUR  wEISmTj  ARE  RETURNED  IN  SER 


COMPILER  I OISS: II 

SUBROUTINE  SENFLTIM.JA.jR.Xt  F.SERI 
WR11EIs.4SSSIm.JA i jR , jc .f 

44(4  fCRNAlt'  SENFLIi  M:«,I4.w  O:',!!,'  JR:*.Ij, 

w  «  jCi‘,IJ,»  Fs* , 612.4> 

DIMENSION  SERI  1 1 

H:j.H  1st  24 S 

ff:M 

DO  100  III ,N 
F1K 

L-lf 1-l.l/FM 

60  101  NO, SO.  40, 41. 411.  JC 

(I  If  IU-. OOOOOOOl  140.40.63 

(1  If  UC-4)44  .44.4S 
(4  <er in-  siNipiwut  yipjwut 
(0  10  TO 

(S  SCR|ImSIWIPl*ul/lPIwul»w,2 
CO  10  10 
40  SERI1K1.-U 
CO  10  TO 

(0  ScRIimSwIi.wCOSIPTWUtl 
CC  10  To 

*0  If  IABS(UI*.S»S1,S2,S2 
SI  jER  1 1 1  =  l. -4.WUW*2«4 , WARSI UW wj I 
CO  10  TO 

<2  SER I I t:2.w |1 .-SRSIUI Iwwl 

TO  IMF-. 000000011100, 100, T1 

11  If  IARSIP1  wlFI-1.  IWFI  -.000011100. 100.  TI 

Tl  Sf  Rll  >:$£R  1 1 1  »S  iNIPI  41  FI  - 1.  |wF  1/1  PI  sir  1-1 .  IWF I 


1 


continue 

iruou  ».  uc.ioi 

IW:4*itt>ll 

•  «  104  |!li» 

II 

CO  lit  1=1 

irui-iitttfUMtii 
i  scatii=-seRii> 
inti  ni».UM» 
i  umiku *u*n  • 

i  continue 
g  -*o=J»*2 

to  ion  sOtiso*i*o#iooitJ*a 

a  continue 
ft  tun 
o  tlM/l 
Cc  iss  1=1 (BB2 
i:SCRIlt 
II=n-l*l 
uatlusentiii 
s  utiinn 

cetwttN 

o  *ns2«n-t 

on  net  otooooi  1 1 «*c*t  iiti=iiWti 
0  roBBOT 41* ,012. Of 
co  io*  l=n,t.-i 
ll=l**-l 

s  hr  ill  i=se Blit 
r«  t  =n- 1 

cc  1*1  l=i*nni 
11=00-1*1 
it  set  4tl=lt*IIIt 
citutn 

ICO 


*10t,S  1f*UC1*P**0.B£H*£ 


11*U0 ?*1R60 

1 

2 
1 
« 

* 

t 

* 

10 

It 

12 

tl 

1* 

IS 

n 

u 

i* 

i* 

20 

21 

22 

11 

2* 

25 

2* 

21 

2* 

2* 

10 

11 

12 

11 


m.ccnivE 

CC*PH.tRtOU6  =  }l 

subboui me  oe*i»eio»iniT| 

OlnCNllon  01201  «*n01  201 
common  /UM»T/nino, 1*01201 
ccBPon/s*Linc/KSfL 

c* i*  cr.iSRl* /i..l/  .  .  ....  ..... 

L001CCL  IOB»,IOM, 1010, IQS*  tlOOT.lOSlB 
C 

1=0121 
1=0111 
C  =01*1 

CO  10  1 21 ,221 tlSBlB 
C 

21  *=0 

IS* 11=2 
CC  I  I:l>2 C 

JMlNO'll.tO.OI  60  TO  1 

It»Sl=I 

»:*«! 

«*0t* 1=1 

ioe«=i  i.eo.s.  or.  ioo*» 

ics«  =  ii.eo*v>OR*ios*i 

14)  1:1  l.CO.IO.OR.lOOtl 

101t:ll.eO.II(OB.I.eO.I2.0*.101T» 

lCF0  =  ll.t0.l2.0R.I01DI 

l0S16  =  ll.tQ.6.0"(l*e0(  1(0*(W*161 
1  continue 
c 

MS*:* 

If  illlSUMIlSi:) 
lflltfSt.CQ.Sl  *Sw:2 
1M|0SI6I*SU=1 


)« 


IT  tIOSIt.ANO.IIlAST.tC.TI  •  KSM  -$ 

c 

22  SO  IS  iSSiSIiMiIMIIiMII 

c 

ll  (0  10  ISII,lli|,ttU 
SI1  «u  OCNOPNIP.I.C.SI 
(0  10  SIS 

s  12  call  KiiitM.CiSi 
si)  msi-s 
c 

SO  10  IS0,S0.S2,S1.S2I,KSM 

c 

*2  cm  siMio.iiSiiiiMin 

CASI-SIOT 

0I)I=SIC 

c 

so  lo  CSO.Sg.SO.Ss.sSI.OSM 
c 

S)  lFIloSVICAlL  SOUNOIt  «S.P«  SO  I 
0I«>:S« 

IftlOPIlCAll  POUNPIP.T.S.TPOTI 
0lll|:TPOT 

IF  I I0P0 1C  ALL  SIONAtP.TPOT  ,S  .SI  SPOT,**  I 
Cllil:IISFol 

IFlIOOSICA  U  ORNO  ALII  ,P.S  .0.10  IT  .CM  .00  ION  I  ,110  V  I 
CIIIUA 

IF C 1001 ICALL  TOUOCMP.SlC.lNll.OCPI 

citeisocp 

c 

SO  C0NI1NUC 

c 

CO  2  1:1, NINO 
K-KND  f 1  I 

2  S 1 1 1-0  IN  I 
C 

KUW 


c 


[NO 


•P01, 


PF*U07Rpt  SO  .NOCCC 


PF*UO TSPOSOC 11*  AOCCC 


1  CONPILCOIOIAO:)) 

2  C 

)  SUOOOUTINC  NOCCC  IP.T.C.CF | 

•  c 

S  C  NOCCC  CORRECTS  FOR  UNPCRATURC  AND  PRCSSUOC  CFF[CTS  ON  CCLL 

S  C  A  AO  APPilCS  CAiTORATlQN  CORRCCIION,  CF. 

T  C 

(  DATA  ALPMA/t.SoC-t/.FK/tT.at/ 

*  C 

10  C  :CF*C*|I  .-AIPHA*|T-1S.  MIPS)  .RSOSOI/FRl 

11  C 

12  RCKON 

I)  CNO 


OPRT.S  PFPUOTRpSSO.OCNOPM 
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WlMHMWlfl.iUlN 

1  CORflLCR(OI*f :ll 

2  C 

s  uutuiiK  KMmirttiiitbiMfi 


« 

t 

* 

t 

i 

« 

10 

It 

12 

11 

If 

If 

If 

IT 

M 

If 

20 

21 

22 

21 

2* 

» 

1* 

21 

20 

2f 

10 

It 

12 

11 

If 

IS 

If 

11 


C  IMS  Routine  out  INC  01  until  11  fit  I  •  «Ht  OtOfttfUKItlth 
c  tot* lions  10  COURUTC  SUIHItf  fROU  RRCSfUOC.  TtnfC**lUR£, 
c  »»o  conouciifiir. 

C  flsfRCSSURC  in  OCCtOROf 

C  IlMCnKRMURC  in  ocoocct  atsiut 

c  ctsconouciiom  in  uu mo/cm 

c  RfRinsconouciiviii  rotio 

c  confute  conouctifttt  n*tio 

RfRlRsCl/fl.tO* 

c  connect  nftin  to  2tno  fnciflioc 

r:i(n.i*ft-tS*ti-s.fofSC-ioi*fi«i.M>oif[-fi*fi/ 

1  <l].im-fftl*1.010M*2lf1|«t.» 

Rsi:ofoinm.«n 

c 

c  connect  nst  r on  tiufcn«tune  errcctf  «t  lifii  nnene 

c  mu:c  1S.1.0/C  1S,IS,0 

c 

R Instill  .ISOfOC-tftl-t.ZMRjC -lift  1*1. Iioooe-f »*ll*2.00S2fft-2l 
I  *ll*.tltSS2f 
RS.-RSWRtn 
C 

c  cntcutm  sntinitf 

nsfsllil-l  .S2SII*nS*S»t002* |*0S”lO»OI0Ot*fOf*12.10002lfOS*20«0S0t  1 
|«RS-.QOm 

R  SRSn  I=RS«|RS<  I.  If  Ilf. f  2t -2-. **C-lfll-f.£-J**S>*!l 
1  «ll.2St~f-2.*I-tf1l»*t  11 

iff  t  :R$f  «n$nsni 
c 

rc  runn 

tnu 


ofntvs  fmuatffCRO.fssiR 


fr«uotff*fioi 

1 

2 
1 
« 

S 

• 

1 

0 

t 

10 

11 

12 
11 
If 
1  • 

I* 

11 

IR 

It 

20 

21 

22 

21 

If 

21 

21 

21 

It 

If 


II. MSI* 

C0UfjLCRI0!tCO}l 
sutnouimr  fSSTtif  .t.c.SAti 

c  mis  Routine  uses  me  inicucu  soLimit  sc*tt  iiftoi 
c  e  out  lions  to  confuic  sttiuitt  rttn  tRcssuRC, 
c  Mute Rtiuot  *no  cououcmitt  . 
c  ncce  jounnit  u*n  ittoii 
c 

C  f  I  :RRt  SSURC  1001 

c  tisienrcRRtunc  loet-ci 

c  cisconouciivitt  iNUMO/cm 

c  Rfsconauctivitt  nttio. 

c 

C  ccnfiiic  conouctlfllt  MtiO 
c 

Rf-C/f 2.tif 

c 

»tfHii2.oioe-s*f-*.iioe-io^*f2.»i.**oe-is*f«*i.» 
t  /ll«S.f2St>2ft«f.fkfC-f*l**l.* 
c  f .2ise-i*nf'i.iotc>i*nf*ti 
c 

c  iLfMfin ,i,f»  is  tne  rn*ctionn«L  inRCnsc  in  cono  owe  to  rrcss. 
c 

SnlC0st.l**0fie-l«2.005**t-2ft«l.l0f21tt-fftf*2. 
e  -*.t*»ot -i*tf*i. *1.00111 -t#t*ff . 
c 

c  mis  is  t  cns.t.ai/ciis.is.oi 
c 

RRtsRf /ISUtCOfl l«**LfHII 


C 


10  S41:|0.00i0-0.IST2*RRT*40.S*IS.mi4RRT 

si  (  «««it-is. 1/11. *10.01*2*1  t-is. tin* 

J]  1  IO.OOOS-O.OOS**ial**0.|-O.OQ***Ml 

]4  t 

IS  (  -0.0l444MT442.su 

I*  t 

IT  OC TORN 

JS  LNO 


•M1.S  PF*U0T*pSSO.Sl6N4 


IT4PRS0I II  .SIRS  t 

1  Sut  ROUTINE  tlMtll.TitilllliUlttM 

2  lRI9i-l-J.TR 

j  S16SlS:t-tNlTi*TNSTi*l2iS.*TII/IS01.ST*IT**?.2SII 

4  SlSOrS*lS*l  .*T*Ti*14C-S*S-. 4S24T* 14(-I  |..R14iToSi> 

5  1  -.9I44SRR  J2C-J 

i  «i:i4l.4Titfr-2-Ta(.filitC**-T4.10MK-f  II 

T  IUI4|.1|0 101-4-1*1 .11 S4E  -i-t*  .  IRRIC  -Til 

•  Slit:  SI6SIS.|SisO«.|J2«|4||.-iT*iT*lSliO-.I  I24II 

4  C  INIS  IS  SISNi-t 


10 

SIS Ji-SlSO-21. 

11 

IIMR  t./l  l.*l.iJE-S*PI 

12 

l2:22T.«T4t2i.  JI-T4I  .SS1-T4J04II 

II 

1 J:i.C-«*P*llOS.S*t4(T.SO-l*. ISI-I.SE-44PII 

14 

14:.14SIS2i4||4T.  J-T4I2.T2-0.044TI 

|S 

1  -1  .E-44R4 112  .4-T4|.4T -.02*11 1 1 

14 

l&:|.C-2*<SlS2i*Sl«2il*l4,S-.t*T-P*II.RC-*-.*C  -S4TII 

IT 

«NU: 1.1-4*11 1-12. Il-14.il  1 

It 

UF*:tl.-INU*Pt/ISt«T4.001*  |.l 

IT 

RHOSIPM./ilFi 

20 

SICS1F:«RNOSTp-I.I41oO0. 

21 

C 

TMS  IS  SUM -SIP 

22 

RlluRN 

21 

ISO 

NPRI.S  PF4U074P  1 00. SC  UNO 


1  SUBROUTINE  SOUNDIT.SU. PRCS4TCLI 

2  C  Iidlsto  Itfssuit  tONTtaSlON  JWO 

1  C  REFERENCE  -  NILSON,  H.0..I940.  EOuMION  FOR  Tilt  SKEO  OF 

4  c  stui  in  sra  m*ter,  jour,  icoust.  sec.  **er.(  jfiiii.iist. 
s  c 

4  C  ms:  FKSSURC  IN  DECItiRS  FROM  S|i  SURF4CE  . 

T  C  F  :  TRIM.  PRESSURE  IN  R0/CN2  RRSOLUTE  . 

■  C  SU  :  SUINTTT  IN  PiRTS  FE|  1000  . 

f  C  1  :  TENFERi  TORE  IN  OEfREES  CELSIUS  . 

0  C  S0U40  :  SOUNO  VELOC ITT  IN  «TERS  PER  SECOIS  . 

1  C 

2  F  :  IFRES  •  I0.II2SI  4  O.IOITTU 

J  S  :  ML  -IS. 

1 4  C 

5  41  :  T*  |4«ST21*T4|-4  .4SI2E-2*T  41 -2  .404  SE  *4*?  .TislE-i*T  III 

4  4}  :  4*11. JTTTT«I.*T202C-J*SI 

T  4F  :  F4(o.|i02T2«F4||.02»iE-t*F4II.S»liE>T-I.J»OJE-124rill 

It  t 

T  9S1F  :  S4't4ft.|244t-2*T.TniE-T41l4P4IT.T0!»E*S-|.2T4lE-T4P« 

!0  1 141  ]. I  SIDE  -i* I .STTOE  -T*f  I  H  * 


29 


21 

2P«t!k<-l.tkOTC-N 

22 

1FNRNI 1*1*2 .S2NNC 

21 

C 

2N 

SOUNO  :  IN  NO .IN 

25 

c 

2k 

ItCL  I  SOUND 

21 

ACT  URN 

20 

CkO 

RPRT.S 

PFkUOTNP  ISO. iRuOCP 

rrwoKNcoi  ii.iracir 


1 

2 
1 
N 
5 
k 
T 
• 
« 

10 

11 

12 

U 

14 

15 
Ik 
IT 
It 
I* 
20 
21 
22 
21 
2» 


tOI#IUItoUt:ll 

SlRROUTtNC  TRUDCP|PRCSS,S1RSTP,INIT,TRU0P| 
CONhOm/MOoCI/IOoCWIII 

connon/rdocf /rose* i I i 
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